LPTHE-98-24 



Structure of Matrix Elements in Quantum Toda Chain. 



F.A. Smirnov Q [] 



Laboratoire de Physique Theorique et Hautes Energies q 
Universite Pierre et Marie Curie, Tour 16 1^'" etage, 4 place Jussieu 
75252 Paris Cedex 05 France 



Abstract. We consider the quantum Toda chain using the method of separation of variables. We show 
that the matrix elements of operators in the model are written in terms of finite number of "deformed Abelian 
integrals" . The properties of these integrals are discussed. We explain that these properties are necessary in 
order to provide the correct number of independent operators. The comparison with the classical theory is 
done. 



OMembre du CNRS 

^On leave from Steklov Mathematical Institute, Fontanka 27, St. Petersburg, 191011, Russia 
^Laboratoire associe au CNRS. 



1 



1 Introduction 



As it became clear recently g, ^ there is a close connection between the formulae for the matrix elements 
in integrable field theory (form factors) Gj and the method of separation of variables developed by Sklyanin 



The form factors are typically given by certain integrals. These kind of formulae can be interpreted as 
follows. Consider an integrable model which allows the separation of variables. The separated variables nat- 
urally split into two equal parts: one of them can be considered as "coordinates" and another as "momenta" 
(of course they have nothing to do with original canonical variables in which the model is formulated) . The 
formulae for the form factors are understood as matrix elements written in "coordinate" representation, i.e. 
in terms of integrals with respect to the "coordinates" . 

Another observation made in and used intensively in |^ is that the integrals in the formulae for 
the form factors in models with sl{2) Lie-Poisson symmetry (Sine-Gordon, for example) can be considered 
as deformations of hyper- elliptic integrals. This fact must be also related to the method of separation of 
variables because the "coordinates" describe classically a divisor on the spectral hyper- elliptic curve. The 
important conclusion made in |^ is that these deformed hyper-elliptic integrals must have similar properties 
to the usual hyper- elliptic integral in order that the correct number of equations of motion exists in the 
quantum case. 

In paper we performed the quasi-classical analysis of the matrix elements in CFT in finite volume. This 
is a much more complicated case than the case of infinite volume. The method of separation of variables 
seems to give the only possible approach to the calculation of form factors. The main difficulty of the 
problem in finite volume is due to the fact that the separation of variables leads to the Baxter equations 
whose solutions describe the wave-functions in "coordinate" representation. So, one must consider integrals 
of solutions over Baxter equations. 

In the present paper we consider a much simpler model which nevertheless exhibits difficulties similar 
to those of integrable field theory in finite volume. This is the periodical Toda chain. Historically this is 
the first model to which the method of separation of variables was applied In this case the problem 
of describing the spectrum leads to Baxter equations with non-trivial solutions in entire functions. The 
matrix elements are given by integrals over these solutions. We show that these integrals can be considered 
as deformed hyper-elliptic integrals allowing deformation of all the important properties of hyper-elliptic 
integrals. Similarly to Q these properties are needed for the correct counting of operators, they are actually 
equivalent to the equations of motion. 

2 Classical Toda chain 

The periodical Toda chain is described by the Hamiltonian: 



i 




(1) 



where Pj, qj are canonical variables, Qn+i = Qi- 

The exact solution is due to existence of Lax representation. Consider the L-operator 




and the monodromy matrix 




Obviously, detM(A) = 1. The monodromy matrix satisfies Sklyanin's Poisson brackets: 



{M(A) ? M(m)} = [r(A - fi), M{X) ® M(m)] 



(2) 



where 



r(A) 



P 

a" 
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P is the permutation. The coefhcients of r(A) = trM (A) are in involution: 



{r(A),T(M)} = 



Moreover, 



T(A) = A" - PA" 



-P^ -H] A""2 



where P — ^Pj is the total momentum and H is the Hamiltonian (|^). Thus T(A) generates n integrals of 
motion in involution providing complete integrability of the system. 

From here on we can forget about the Toda chain saying that we consider an orbit of Lie-Poisson group 
i.e. the polynomial matrix M(A) with det Af (A) = 1 satisfying the Poisson brackets (||) (determinant is in 
the center of these Poisson brackets) and characterized by certain reality conditions which we shall discuss 
later. 

Let us consider the elements of M (A) in some more details. We introduce the notations 



A{X) = A" + A''' -ai 



+ A^ 



B{\) =b{X 
C{X) = A"-ic2 
D{X) ^ A"-2^2 

The variables b and ai have the Poisson brackets 



n-l 
n— 1 I \n— 2 



&i + --- + 6„_i) 
+ c„+i 



{ai,b} = b 

and Poisson commute with the rest of variables. In terms of Toda chain ai = P and b = e"^" describe the 
motion of the center of mass. Our nearest concern is the algebra of observables A. We define this algebra as 
the one generated by all the monomials of finite degree of the variables aj , bj , Cj , dj and b. It is important 
that the polynomial structure of M(A) introduces grading of A. Namely, we can prescribe the degree i to 
every of elements a^, bi, Ci, di and degree to b. The degrees of the leading coefficients of C(A) and ^(A) 
are chosen in order that the coefficients of the determinant 



detM(A) = A2"-V2 + 



/2 



are homogeneous. The variable & is a kind of zero-mode, it is of minor dynamical value. The algebra A 
contains a subalgebra of polynomial functions of Oi, &i, di and Ci = bci. So, this subalgebra does not have 
6 as a separate generator, the change in definition of Ci is needed in order that the Poisson brackets are 
closed for A{X), -D(A), B{X) = n(-^ - 7i) and C{X) = bC(X). We shall deal only with this subalgebra. 

The algebra Ao as a vector space splits into direct sum of subspaces of different degrees. Let us denote 
by 6{n) the dimension of the subspace of the degree n. The generating function of S{n) (character) is given 

by 



Xiq) 



oo ^ 



1 



[1] 



n=0 



[1] [2n]\ 



where [n] = 1 — g", [n]\ = [1][2] ■•■[«]. The first four multipliers come from monomials of of aj,bj,Cj,dj 
respectively, the last multiplier comes from factorization by the condition detM(A) = 1. 
Notice that 



X(9) 



where we introduced the q-binomial coefficients 



1 




'2n - 1' 




'In - \ 




[n]\ [n- 


^( 




- q 


n - 2 


) 



(3) 



m ! n — TO 



Later we shall provide an interesting interpretation of this formula. 
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Let us return to more traditional consideration of classical Toda chain. We do not give a complete list 
of references which can be found in |5|, the important for us fact concerning the classical system is that it 
allows the separation of variables 9|. Consider zeros of the polynomial B{X): 

and the variables Aj = D{jj). Notice that Aj = A(7j) where A(A) is the eigenvalue of A/(A). The variables 
7j, logAj are canonocally conjugated which can be shown following [|j using (||): 

{7j,logAj} = (5,, J 

From det Af (A) = 1 it follows that A{'jj) = Aj^. One can reconstruct the matrix M{\) from 71, ■ • • ,7„_i, 
Ai, • • • , A„_i, ai and b. The simplectic form is written as 

n-l 

u) = d log Aj A d-fj + d log b Adai 

The 1-form a (uj = da) is 

n-l 

a — log Ajdjj + log b dai 

Let us take other coordinates on the phase space, namely, 71, • • • , 7n_i, t2, - ■ ■ ,tn (defined by T'(A) — 
A" + X^-^ti + y-'^t2 H h <„), h = ai and b. From 



A, = A(7,)-i(^r(7,) + .yr(T^ 



one easily finds the expression for the symplectic form in these variables 



n — 1 n 'y^^^^' 

Lo = \] / dtk A d'jj + d log b A dai 



where P(A) = T'^(A) — 4. Thus the equations of motion take the form 

A - 7fc 



{T(A),6} = A"-i6 (4) 
certainly, only the first n — 1 equations are really interesting. They are linearized by the Abel transformation: 

n-l 



A-'-i 



where aj are the first kind abelian differentials on the spectral curve /i^ = ^'(A): 

dX 



We associate the "times" ri, ■ • ■ , Tn-i with t2, - ■ ■ ,tn: 

A.F^d,F^{t,+,,F} 

n-l 

The evolution of ^ J^*" crj with respect to times is linear. 

k=l 
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The above considerations apply to any orbit of the Lie-poisson group. We want now to consider specific 
reaUty conditions which correspond to Toda chain. It can be shown |^ that the conditions in question 
are: 

1. The polynomial T{X) of degree n has n real zeros. Moreover its local maxima are not below 2 and its 
local minima are not above —2. So, all the zeros of the polynomial P(A) are also real, they are denoted by 
Ai < A2 < • • • < A2,i. 

2. The polynomial -B(A) has real zeros 71, • • • , 7„_i which belong to the "forbidden zones" : X2k < Ik < A2fc+i. 

The equations of motion (^) preserve these conditions. The hyper-elliptic Riemann surface = P(A) 
has 2n branch pints (A^). Its genus equals n — 1. We present the surface as two complex planes with the 
cuts along (—00, Ai], [A2, A3], ... ,[A2n,oo) identifying the banks of the cuts on two sheets in usual way. The 
canonocal a-cycles Uj are taken as ones encircling the cuts [A2j, A2j+i] for j = 1, ■ • ■ , n — 1. Topologically 
the points jj move along the cycles aj. 

Define the normalized holomorphic differentials 



such that 



Then 



j — Ajk o-fc 



fc=i 



are real angles on the Jacobi variety, and the dynamics describes linear motion along this real torus. One can 
invert the Abel transformation expressing the symmetric functions of 71, ■ • ■ , 7n_i (recall that they coincide 
with 61, • ■ • , hn-i) as functions on the Jacobi variety (functions of 0's) using the Riemann theta-function but 
we shall not need explicit formulae. The angles 9 and the times t are related linearly: 

n-1 

1=1 

so, the using the theta-function formulae mentioned above one can resolve the equations of motion expressing 
bj as bj = bj{Ti, ■ ■ ■ ,r„_i). 

From the point of view of algebraic geometry the monodromy matrix M{\) gives an affine model of 
hyper-elliptic Jacobian, and the functions 6j(r) are generalized Weierstrass functions In the case of 

genus one (n = 2) the function 7(ti) — 6i(ti) is the usual Weierstrass function which satisfies the second 
order differential equation 

9h = ~ni) (5) 

One of results of our further analysis will be in finding certain second order partial differential equations for 
generalized Weierstrass functions which can be thought about as generalizations of (||) . 

Let us consider the ring of generalized Weierstrass functions with coefficients in fi, • ■ • , i.e. the ring 
of polynomials 

F{ti, ■ ■ ■ tn,bi, ■ ■ ■ , bn-l) 

Consider further all possible derivatives of these polynomials with respect to r^: 

dt---dt'riFihr--tn,bi,---,bn-i) (6) 

The equations we are looking for correspond to all possible linear combinations of the functions (^ which 
vanish due to equations of motion. To understand the origin of these equations we have to return to our 
mechanical considerations. 

Mechanically one understands the derivatives di as action of hamiltonian vector-fields. Using the Poisson 
brackets (||) one can express (^) as a function of ai, • ■ • , a„, 61, • • • , bn-i, £2, • • • ,Cn+i, ^2, • • • , rfn i-e. as an 
element of the algebra ^0 ■ We put forward the following 

Conjecture 1. Every element of Ao can be presented as linear combination of the expressions (|^). 
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We were not able to find a complete proof of this statement, however, the consideration of examples supports 
it. Further indirect support of this conjecture will be provided by the calculation of characters given below. 

Assuming that the conjecture is true one realizes that the way of presenting an element of as a linear 
combination of the expressions may be not unique. Indeed, let us calculate the character of the the space 
span by (^. We prescribe the degree i to di which is consistent with the Poisson brackets (||). Obviously, 
the character is 

7Z TTTIZTTTZ TTT > x{q) 



where xil) is the character (||). So, there must be linear dependence between the functions (||) which is 
responsible for differential equations on the generalized Weierstrass functions. Moreover, there is a criterium 
which allows to judge whether the set of equations is complete. Indeed, to show the completeness of the 
equations one has, obviously, to prove that taking them into account leads to correct character (^). 

Let us find the equations in question. To this end we shall use the Fourier transform. Consider a function 
F{ti, ■ ■ ■ tn,bi, ■ ■ ■ , The variables tj are the integrals of motion (and the moduli of the Riemann surface) 

and the variables bj are the functions on the Jacobi variety due to the equations of motion. Hence 

F(ti,---t„,5i(T),---,6„_i(r)) = 

= J2 j de[--- j de'^_,F{h,---tMo')r--.bn-i{o')y'''''^ 

where dj = AjiT^-i- Let us undo the Abel transformation inside the integrals: 
•••*„, 61 (t),---,6„_i(t)) = 

xi?(ti,...i„,7i,...,7„_i)ne^*'=^^^-^ (7) 



where F{t,^) = 6(7)) (recall that bj are elementary symmetric functions of 7's), for any 
k = {ki, • • • , fc„_i} we define 



Ml) = Jt 
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Deriving (^ we needed to take into account the Jacobian of the Abel transformation. Later we shall relate 
the integrals in (^ to the quasi-classical limit of the matrix elements in quantum Toda chain. The equations 
of motion correspond to vanishing of all the integrals in (|^). Let us explain the possible reasons for these 
integrals to vanish. 

Consider first the term in (^ with fc = which is nothing but the average of F over the Jacobi variety 
(the trajectory): 

There are two reasons for this integral to vanish. First one is obvious, it is due to existance of exact forms. 
With arbitrary polynomial ^(7) associate the polynomial 

We marked explicitly the dependence on the moduli (integrals of motion) t = {ti, ■ ■ ■ ,tn} which enters 
through P(7). 

There is an obvious 
Proposition 1. The following integral vanishes 

^ A(i)(7) = (9) 
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for any polynomial L and any closed cycle c. 
Hence the integal m) vanishes if 



n-l 



^(*'^) = E TtT^ T A(L)(7.) G(7i,---7.---,7n-i) 

^ n(7. -7j) 

for any polynomial L and any symmetric polynomial of n— 2 variables G (both of them can be also polynomials 
of parameters t) . This property means in particular that by adding exact forms one can reduce the degree 
of the polynomial F in every "fj up to n. 

The second reason for the integral to vanish is due to Riemann bilinear identity. Consider the 
anti-symmetric polynomial of two variables 

Ct(7i7 72) = ^t(7i:72) - -Rt(72,7i) (10) 

where 

i?t(7i,72) - V^(7i)/- (^—VpM 
«7i V 71 - 72 

For any two cycles on the Riemann surface one has 

C't(7l;72) = Cl O C2 



where o means the intersection number. Since the cycles Oj do not intersect one has the following 
Proposition 2. For any two a-cycles aj and Ok the following integral vanishes: 

^ ^ Ci(7i,72) = 



/P(70 /P(70 
Hence the integral (^) vanishes if 

7) = 7 ^ T-r / V( V Ct{-ii, 7j) G(7i, • ■ • Ti ■ • ■ 7j • ■ • , 7„_i) 

fr^ (7* - 7j) n Kli - liKlj - li) 

Let us consider now the case of arbitrary k = {ki, ■ • ■ , A:„_i}. Introduce the polynomials Sty. 

i kj uj^M — —^=^21^7 

Integrating by parts one gets the following three simple propositions 
Proposition I'.For any polynomial L define the polynomial 

7 

DtML)h) = A(i)(7) - StAl) J L{^')StAl')di 



then the following integral vanishes for any a-cycle: 



^!2_ e'^'°(^) Afe(L)(7) = 



Proposition 2^. Define the polynomial 

71 72 

L't,fe(7i)72) = Cf(7i,72) - ifc(7i) / a7 + '^'fc(72) / oj 

J 7-72 J 7-71 
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then for any two a- cycles Oj and Ok the following integral vanishes: 



d'fi 



VWii)J /P(72) 



There is one more identity which is trivial in the case = 0: 
Proposition 3'. For any a-cycle the following integral vanishes: 



From these three Propositions one finds the foUowing partial dilFerential equations on the symmetric 
functions of 7: 

1. For any polynomial one variable L and any symmetric polynomial of n — 2 variables G the equation 
holds: 



E 



1 



A(i)(7.) G(7i,---7;---,7n-i) 



n-1 



- ^ didr, 



/,m— 1 



E 



1 



-i n(7.-7,) 



7,; 



"d7 G'(7i,---7j---,7„_i) 











(11) 



2. For any symmetric polynomial of n — 3 variables G the equation holds 

1 

Tl-l 



C(G) ^ 5: 



(7i -7j) n (7^ -7i)(7j -7i) 



Ct(7i, 7i) G(7i, • • • 7,; • • • 7j • • • , 7„_i) 



5] a,a„ 



,m— 1 



E 



1 



(7« -7j) n (7i -7i)(7j -7i) 



G(7i, . . .fi . . .fj . . ■ ,7„_i) 



X 7 



n-l-l 



,,n— 1 — rn „,n— 1 — m 



7j 



7-7j 



-^7 — 7! 



i 



7-7» 



c?7 



= 



(12) 



3. For any symmetric polynomial of n — 2 variables G the equation holds: 



ri-l 



E 



1 



7r'"'G(7i,---7^---,7«-i) 



~[ 0(7.-7,) 



= 



(13) 



One must add to the equations (|11|), (|12D, (|13|) their trivial consequences i.e. the equations obtained from 
them by applying arbitrary number of di. We claim that in this way all the equations of motion can be 
described. 

Let us illustrate this point considering the simple example n = 2. In that case genus is equal to 1 and we 
have only one variable 7. The equations (p^,(p^ are trivial, so, we are left with (pi]). Let us take L(7) = 7P, 
p = 0, 1, • • •. Then (|Tl|) turns into 



Pi"-' P{l) + \y ^P{i) = dl f ^7' 
2 d-y VP + 1 



(14) 



we have only one time n in that case. As it has been said earlier the function 7 is the Weierstras "P-function. 
The equation (|lj) coincides with the usual equation (^) when p = 0. For other p we get the equations on 
degrees of 7 which can be verified for the Weierstrass function. Recall that we were considering the space of 
functions of the kind 



d^Fiti,t2n) 



(15) 
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Let us calculate the character of this space modulo the equations and their trivial consequences (those 
obtained by application of di). Obviously, using (Q) we can express any function of the kind (^5|) as linear 
combination of 



So, the character is 



Q Q 



1 

m V ' [1] ' [1], 

which coincides with (||). 

Let us consider the general case. We have the space of symmetric functins of 71, ■ • ■ , 7^-1 with coefficients 
in ^1, • • • , t„ on which the derivatives 9i, • • • , 9„-_i act. We can consider the derivatives also as coefficients 
identifying this space with -ff„_i which is the space of symmetric polynomials of 71, • • • , 7„_i with coefficients 
inti, ■ ■ ■ ,tn and di, ■ ■ ■ , dn-i (recall that di and tj commute). Substracting the exact forms ( pi] ) we finish with 
the space which is the subspace of defined by the condition that the degree of the polynomials in 

every jj does not exceed n. We define the spaces Hj {j < n ~ 1) of symmetric polynomials of j variables 7i 
whose degree in every variable does not exceed 2n — j — 1 with coefficients in 1 1 , • • • , i„ and (?i , • • ■ , (9„_ 1 . The 
action of the operators C and Q defined by (|l^),(^3|) can be obviously extended to the action from Hn-3 to 
Hn-i and from Hn-2 to respectively. It is also clear that the images of respectively i?„_3 and i?„-2 

belong to i?n-i- One can easily generalize the definitions of C and Q allowing them to act respectively from 
Hj-2 to Hj and from Hj^i to Hj. For these operators one has 



[C, Q] - 0, 



Q' = 0, 



Kerl 



Hi -2 







(16) 



Now noticing that deg(C) = 2 and deg(Q) = 1 one evaluates the character: 



2n - 1 
n ~ I 

n — 3 
1 



2n- 1 
n-2 
'2n - 1 
n - 4 
2n - 1' 
n — I 



2n - 1 
n — 3 

.2 



q 



2n 

n-2 



2n- 1 
n — 5 
1" 



-■■■ + i-qr-' 



2n-l 


'2n - 1 




[n]\ [n- 1]! 

which coincides with the character (^) . There is obvious similarity of what we have done with the paper 



3 Quantum Toda chain 

Following Sklyanin |^ we shall use the same notations for the quantum analogues of classical objects that 
have been used in the classical case. The Hamiltonian of the quantum periodical Toda chain is given by (Q) 
with pj , qj being the canonical operators 

[pi,qj] = ifi^i,] 

Consider the same definition of L-operator and monodromy matrix as in classics. The monodromy matrix 
satisfies the relations 

R{\ - Ai)(M(A) M{ii)) = {Mill) ® I){M{\) ® I)R{\ - n) (17) 
with i?(A) being the quantum R-matrix: 

i?(A) ^X-ihP 

The trace of the monodromy matrix provides n commutative integrals of motion. The center of the algebra 
is created by the quantum determinant 

A{X)D{X + in) - B{X)C{X + in) ^1 (18) 

The idea of using the separated variables in quantum case goes back to . It was developed as universal 
method by Sklyanin. Let us briefly review the method of separation of variables following From the 
relations ( p7| ) one finds, in particular, that 

[B{X),B{fi)]=Q 
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So, presenting the operator B{X) in the form 
one gets a family of commuting operators: 

As it was in the classical case the operators b and ai = P commute with everything except between them- 
selves: 

[ai, b] = ihh 

The idea of the method of separation of variables is in considering the model in b, 7-representation. The 
canonically conjugated operator to b exists already: this is ai. The canonically conjugated operators for 7^ 
are constructed as follows. Consider the operators 

A{X) = A" + A"-iai + --- + a„ 
D{\) = X"-^d2 + ■■■+(!„ 

then it can be shown that the operators 

Aj =7" + 7""^ai + --- + a„ 

A, =77-2^2 + (19) 

satisfy 

A,A, = 1, [A„7g] =zM,-feA, (20) 

The order of multipliers in ( [l9| ) is important. 

It is possible to reconstruct the operators A{X), B{X),C{X), D{X) using ai,b,jj, Aj. In particular, 

n— 1 ^ . N n— 1 

^W-En(:^ (A. + A-) + n(A-7,)(A + ar + E7.) (21) 

k=i j^k 'J ^ 

The Hilbert space splits into a direct sum of orthogonal subspaces Hp corresponding to different eigenvalues 
p of the zero-mode ai. Let us consider the space Hq, the solution in other subspaces Hp are obtained from 
the one for Hq by simple transformation. In Hq the eigenfunctions of T(A) with the eigenvalue t(X) in 
7-representation can be looked for in the form 

n-l 

( 1 1 7 ) = n ^?(7^) 

i=i 

Applying (|2^) one finds with the following equation for Q{j)- 

th)Q{i)^Q{i + in) + Q{i-in) (22) 

where t{X) is the eigenvalue of T(A) on | i ). In the subspace Hq this eigenvalue must be a polynomial of 
the kind 

t{X) = A" + 0(A"-2) 

The equation ( p2|) is one equation for two unknowns: t and Q, so, at the first glance it is rather useless. 
However, assuming certain analytical properties of (3(7) this single equation actually defines the spectrum. 
Namely, require that the function (5(7) is an entire function of 7 with infinitely many real zeros. Moreover, 
impose the following asymptotic: 

(3(A) - cos ( ^log ( - ) + 7 ) , A -00, 



4 



An 



vr 



Q(A)^e^^"cos( :^log(^-^j +- ) , A ^-00 (23) 
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The normalization of Q is not the same as in For n = 4fc the function Q differs from the function ip from 
1^ by the multipher exp(7rAn/2?i) which is in this case a quasi-constant: does not disturb the equations (p2[). 
If n 7^ 4fc there is essential discrepancy with however, we claim that these are (^3|) which describe correct 
asymptotical behaviour. Our normalization basically coincides with the one accepted in [ p^ . According to 
it provides the only way to have correct quasi-classical limit. We are going to discuss this limit in the 
next section. 

The main conjecture of the paper [|| is that the spectrum of the model is defined by all the solutions to 
the equations ( p^ ) with the analytical properties of t and Q described above and the asymptotic behaviour 
of Q given by the formulae (p3[). This conjecture was proven by Gaudin and Pasquier p^ . 

Now we want to discuss the properties of the matrix elements of the operators. To consider the matrix 
elements we need to know the scalar product in the space of functions of 7j . This scalar product was found 
by Sklyanin |H, let us repeat the essential steps because again there will be a minor difference if n ^ 4fc. 
Consider an operator O which is given by a symmetric function -^(71, ■ ■ ■ ,7,1-1). The wave- functions are 
real, so, the matrix element is given by 

°? °? n-l 

{t\0\t')^ j d7i... j d7„.i J]g(7,)g'(7j)F(7i,-.-,7„_i)u;(7i,...7„_i) 

— OG — 00 J — 1 

where w is certain weight. Requiring that the operator T(A) for real A is self-adjoint and using the formula 
( pl| ) one concludes that 

w(7i, • • ■7„-i) ^n^^* -7i)^«(7i,---7n-i) 

where {5(71, • ■ • 7„_i) is i?i-periodic entire function of its arguments. There are two formal reasons for 
choosing particular w. First, everything under the integral is symmetric with respect to 71, • • •7n-i except 
for the multiplier Wi^j in so it does not make sense to put in w anything that can be killed by anti- 

symmetrization. Second, requiring convergence of the integrals one finds from the asymptotic ( ]2^ ) that w 
can contain exp(— 2777^771/ ?i) with 1 < to < ^ 1 • These two requirements fix w up to anti-symmctrization: 

Tl-l 

^^(7i,---7n^i)= n^*'"''^'""^ (24) 

Again, if n = 4A: this essentially coincides wit h p| , otherwise there is a minor discrepancy. 

Informal reason for taking w in the form ( |24| ) refers to the quasi-classics which will be discussed in the 
next section. 

Similarly to the classical case the algebra of observables ^0 is defined as the algebra generated by all 
the coefficients of ^(A), B{\) = h-^B{\), C{\) = b C(A), D{\). The commutation relations @ and 



the quantum determinant (18) provide sufhciently many relations to show that the quantum algebra of 
observables has the same size as the classical one. To make this statement mathematically rigorous one says 
that their characters coincide. 

Now we formulate an analogue of the Conjecture 1 of the previous section. 
Conjecture 2. Every quantum observable O can be presented in the form 

O^GLitir- ■,tn)F{bi, • • ■ , 6„_i)Gfl(ti, ■•■,*„) 

where Gl, F , Gr are polynomials. 

This conjecture looks more natural than its classical counterpart and explains the mystery of the latter. The 
point is that there is no closed formula for the commutation relations of T(A) and -B(/i) which would allow 
ordering of the operators ti and bj . 

Taking the matrix element between two eigenstates of the Hamiltonians one can essentially neglect the 
polynomials Gl and Gr because acting on the eigenstates they produce the eigenvalues. Hence we are mostly 
interested in the matrix elements of the operators Oq — F{bi, ■ • ■ , 6„_i) which are given by the integrals of 
the king 

°? °? n-l n-l n-l 

( i I 1 1' ) = y d7i • • • y d^n-i n q(7.)q'(7.) n(^» - ^^■) n ^(^i- ■ ■ ■ ' ^-i) n e^^^-^^-) (25) 

-00 -00 J=i i<i J=i J=i 
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where 



F{-fu ■ ■ • ,7n-i) = Fihi-f), • • • , 6„_i(7)) (26) 

The function 

is an anti-symmetric polynomial of 71, • • ■ ,7n-i- Every anti-symmetric polynomial can be presented as a 
linear combination of Schur functions, so, without loss of generality we can replace (|2^) by 

det (F,(7,)) 

for some polynomials of one variable Fi, - ■ ■ ,F„_i. Then the matrix elememt can be presented as a deter- 
minant of (n — 1) X (n — 1)- matrix composed of one-fold integrals: 

det ^ J g(7)Q'(7)^^z(7)e^^(^-"^ j (27) 

We consider the integral 

00 

' Q(7)Q'(7)^(7)e-^^'d7 (28) 



as a deformation of hyper-elliptic integral. The study of the properties of these integrals is very important 
for understanding of matrix elements of operators in the model. 

Since the quantum algebra of observables has the same character as the classical one the equations (^l]) , 
(|T^), ( |l3| ) must have their quantum counterparts. Hence there must be identities for the integrals ( p8| ) from 
which these quantum counterparts follow. These identities can indeed be found, we are going to describe 
them in the same order as it has been done in classical case. 

Let us introduce the operation A which attaches to every function -^(7) the function 

A(F)(7) = F{-f + ih) - F(7 - in) 

The operation A^^ is not always defined, but on polynomials it is. For any polynomial L one can define a 
polynomial F such that ^(7) — A~^{L){j). For uniqueness we require also that A^^(L)(0) ~ 0. Consider 
the integral ( |2^ ) with (5(7), Q'{'~l) satisfying the equations ( p2| ) with the eigenvalues t(7) and ^'(7). For any 
given polynomial L[\) construct the polynomial 

D{L\,,{l) = th)A-\Lt){j) + t'(7)A-i(Lt')(7) - 

-t(7)A-i(Li')(7 - - <'(7)A-i(Lt)(7 - ih) - 

~L{-f)ti-/)t'h) + i(7 + - L{-y - ih) (29) 

Then we have the following analogue of the Propositions 1 and 1': 
Proposition 1" .For any 1 < k < n — 1 the following integral vanishes: 

Q{^)Q'{^)D{L)t,t' (7)e-*^'d7 = (30) 

Using these q-exact forms (|2^) we can always reduce the degree of the polynomial ^(7) under the integral 
( |2^ ) in order that it does not exceed 2n — 2. So, we are left with finite number of different integrals ( [2^ ) 
with F{j) — 7^, j = 0, 1, • ■ • , 2ri — 2. These integrals are subject to further relations. 

In order to define a quantum analogue of the polynomial Ct (71, 72) (pX|). we shall need some preparations. 
Consider the function 

7 — 
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The notation A ^{U{-,6))(j) means that A ^ is appUed to the first argument, i.e. 

A(A-i({/(.,5)))(7)=C/(7,<5) 

The function U' is defined in the same way replacing t hy t' . 
Consider the anti-symmetric polynomial of two variables 

Ct^t'(7l,72) = RtMill,l2) - Rt,t'{l2,ll) 

where i?t,t' (71, 72) is defined as follows: 

i?t,t'(7i,72)=i(7i)A-i(t/(-,72))(7i)+i'(7i)A-i(;7'(-,72))(7i)- 

-i(7i)A-i(C/'(., 72))(7i - iTt) - t'{ji)A-\U{; 72))(7i - 

^ ^ W7i)-i(72))(t'(7i)-i'(72)) 



2 71 - 72 



We have the 

Proposition 2". For any 1 < k,l < n — 1 the following integral vanishes: 

Q(7i)Q'(7i)Q(72)Q'(72)Ct^t'(7i,72)e-*^^'=e-*''='d7irf72 = (31) 



Finally, let us define 

St,Ai)^ti^)^t'i^) 

then we have 

Proposition 3" . For any 1 < k < n — I the following integral vanishes: 

Q(7)Q'(7)^M'(7)e-^^'d7 = (32) 

This is in fact the simplest relation, we consider it last for historical reasons. One obvious consequence 
of this relation is orthogonality of the wave- function. 

Let us say a few words about proof of all these relations. Consider the simplest one (|3^) we have 

Q{i)Q'{i)St,Ai) = Q(7)Q'(7)(t(7) - t'h)) = 

= Q(7 + ih)Q'ij) + Q(7 - ih)Q'i-f) - Q(7)Q'(7 + ^h) - Q{j)Q'{j - ih) 



So, 



Q(7)Q'(7)5M'(7)e-^^'=rf7 = 

00 

Q(7 + ^?i)g'(7)e-^^'=d7 + j Qil - ih)Q' i-f)e-'-^^''dj 

— 00 
CO 

Q{l)Q'{l + iK)e-'-^^^d-i- I Q{-f)Q'i-f-ih)e-^^''d-f 



by shift of contour the first integral cancels the forth one and the second cancels the third one. The shift 
of contour is possible because the function Q behaves asymptotically on the line Im(7) = const essentially 
in the same way as it behaves at the real axis. The relations ( ^0| ) and (|3l| ) are proven similarly, but require 
more sophisticated calculations. 



Using the equations (30), (|31|), ( p2D we can write down equations similar to (|l l\) , (|12[), (|13| ) which would 
guarantee that the number of operators (character of the algebra of observables) is correct. The corresponding 
calculation does not differ from the one presented earlier for the classical case. 
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4 Quasi-classical case 



We had a number of identities described by Propositions 1,1',1"; 2, 2', 2"; 3', 3". They look as reflection of the 
same structure. The goal of this section is to explain the relation between different levels of deformation. 

Consider the quasi-classical quantization of the Toda lattice. We had the following formulae for the 
symplectic form and corresponding 1-form in the coordinates t and 7: 

n-l 

a = ^logA(7j)d7j 
i=i 

uj = da^Y.J2 ^r^ ^^fc ^ dj„ (33) 

i=i fe=2 v^y^j) 

it this section we shall ignore the contribution from the center of mass variables ai and h which is easy to 
find if needed. Using these formulae one immediately writes a quasi-classical proposal for the wave-function 
of the states with eigenvalues t in 7-representation: 



*t(7) = exp (J^J a 

Uh. - 7,)i n (t^)) ' (i r A(7)rf7) (34) 



where fi is defined as follows 

A"-~^{uj) = ^ c?7i A • • • A d7„_i A dt2 A • • • A dt„ 



The relation of the formula (|3j) to the exact quantum formulae that we had before is clear. 

1. The multiplier Y[ ili ~ 7j)^ is a piece of the weight of integration w (|2^). Another one will come from 

the second wave- function in the matrix elements. The wave function without this multiplier will be denoted 
by *t(7)- 

2. The function '^t{l) splits into product of the expressions 

1 \ i f I 



exp J logA(7)d7 



which has to be related to the quasi-classical limit of function Q. One should be careful here because the 
function log A(7) is multi- valued, so, the branches must be defined. Moreover, we could probably need linear 
combination of the wave- functions corresponding to different branches. 
Recall that 

A(7) = I (r(7) + v/T^(7)-4) 

where T{X) is a polynomial of degree n. The polynomial Pij) — T'^il) ~ 4 has 2n real simple roots 
Ai < • • • < Xin- We define the function A(7) on the plane with the cuts along the intervals; /q = (—00, Ai], 
i\ = [A2, A3], /„ = [A2n, cxd). Require that log A(7 ± iO) are real for 7 > A2,i then, obviously, 

log A(7 + iO) + log A(7 - iO) = 0, 76 /„ 

Continuing analytically log A (7) into the plane with the cuts one finds that 



log A (7 ) = log A (7), and log A(7 + iO) + log A(7 - iO) = 27ri(n - j), 7 G ^■ 

The variable 7^ belongs classically to Ij. So, we need the wave function real on entire real axis and not 
containing the factors exp(— ^s) with real s when all the variables are in the classically permited places. 
This requires, first, taking different brances of log A for different 7^ and, second, taking a sum of two wave 
functions with 7^ + zO and 7^ — zO for every 7^. The result is 

*t(7)-ne^^'""^^^^??^(^^-) 
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where 



1 \' /I 



7 



and the branch of log A is defined above, (""pp;;)) " is real positive for A2j-i < 7 < 

To ensure that Qqd"/) is a single- valued function in the plane with the cuts the Bohr-Sommerfeld quan- 
tization condition must hold 

Jj = I logA(7)d7 = nh{2nj + 1) (35) 



the integrals Jj are the classical actions. The cuts of Qqc come from condensation of zeros of the quantum 
Q in the quasi-classical limit. 

Let us consider now the quasi-classical limit ?i — > of the matrix elements. This limit makes sense 
literally if two condition are satisfied: 

1. The quantum numbers are large. In our case it means that the zones [X2j , X2j+i] do not collapse. 

2. The states | t ) and | t' ) are close, i.e. the eigenvalues of Hamiltonians are close: t'j — tj — 0(h). 

Consider the matrix element (|25| ) for such close states. It consists of the integrals 



From the quasi-classical estimation of Q, Q' one concludes that 
Q(7)Q'(7)^^.(7)e^^^-"^^d7 
4 1 _i==cos(^i y\e(logA(7' + *0))d7' + j)cosQ J\e{\ogA'ij' + iO})dj' + F,{j)dj 



We have 



2 cos Q y Re(log A(7' + iQ))di + cos Q ^ Re(log A'(7' + iQ))di + J 
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cos - / Re(log A'(7' -I- iO) - log A(7' -I- iQ))di 



sm 



7 



Re (log A' (7' -I- iO) -I- log A (7' -t- iO))d-i' 



The second term in the RHS can be thrown away because it is rapidly oscillating in the classical limit. The 
variation log A(7) — log A'(7) is estimated as follows. Consider the classical solution with Tc/(A) = <(A) where 
t(A) is one of the eigenvalues. This is the place where using the same notations for classical and quantum 
observables can be misleading, so, we mark explicitly the classical ones. The eigenvalue t'{X) — Tci{\)+5T{\). 
One has (in further calculations we neglect contributions of order o{h)): 

51ogA(7) = log A' (7) - logA(7) - ^^^^^ 



How to find 5T(^)1 The quasi-classical states are subject to the Bohr-Sommerfeld quantization conditions 
(|35|). The quantum numbers rij are quasi-classically large: rij = 0{h^^), but their differences for the close 
states are finite kj = n'j — rij = 0{1). Hence 

6 j log A(7)d7 = / ^^^7 = / ^hTt'^^ = (36) 
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The matrix 

is the inverse for the matrix used in definition of the normalized Abehan differentials . Hence 



dt 



n-l + l 



hk,A 



which means that 



Thus the quasi-classical matrix element for the close states is 



t\0\t')= f^^--- f ^21=^ J](7,-7^.)i^(7i,---,7«-i)n2c°^(<f'=(7.)) 
/ V^(7i) / V^(7«-i) J 



(37) 



recall the notation ^k{l) — P ^j^j- This is the same expression as for the Fourier coefficient in (|^). At the 
first glance there are two disagreements: in we integrate over aj, and we have exp(i<i>fe(7j)) under the 
integral. Actually, these two disagreements compensate each other because aj = {Ij + iO) — {Ij — iO) and 
^kil + iO) + ^kil — iO) = when j £ Ij. Notice that the states | t ) are not normalized. 

Certainly it is not a wonder that we have found the Fourier coefficient as the quasi-classical limit of the 
matrix element. We have performed all the calculations in order to have complete mathematical picture. 
On the other hand from the point of view of physics one can argue as follows. 

Consider the action-angle variables Ji, • • • , Jn-i, Oi, ■ ■ ■ ,9j. The Bohr-Sommerfeld quantization in this 
variables does not give correct quantum result, but still it is correct quasi-classically. Consider the eigenstate 
of Hamiltonians | t ) and the eigenstates of the angles 16*). Quasi-classically one has for the wave- function: 



{t\0) 



J^(^^ 



(38) 



Consider an operator O. This operator can be, at least quasi-classically, ordered in such a way that J's are 
to the left of 6''s. Then the classical shape of the corresponding observable on the solution with given values 
of integral {t) is 

{t\0\9) 



Oc 



h,- - ■ , Sn-i) = lim 



{t\9) 



Insert the complete set of eigenstates into this formula 

{t\0\9} 



lim 



{t\9) 



Y^ {t\o\t'){t'\ 

t' 



) 



{t'\t') {t\9) 



(39) 



Quasi-classically only close states are important for which we have from (38): 

{t'\9) 



{t\9) 



Now it is obvious that ( p9D gives the Fourier transformation of Od (^i , 
(|37|) that quasi-classically 



It is clear from the equation 



det(A) 



{f \t')^{t\t)+0{n), {t\t) 

so, there is complete agreement with the formula (0). 

Let us consider two close states introducing the notation: St^kil) — ^'{l) ~ ^il)- The quantum matrix 
element goes to the classical Fourier coefficient when ?i — > 0. If we consider the classical Fourier coefficient 
with A; = it describes the classical limit of the quantum expectation value {t \0 \ t). So, there is no wonder 
that we have the following sequences: 

A;=0 



A,t'(i)(7) 



Dt,k{L){l) 

Ct,k{ll^l2) 

StAi) 



k=0 
k=0 



Ct(7i,72) 




(40) 



So, we have two levels of deformation of hyper-elliptic differentials. The impression is that the quantum 
deformation is very natural, and that the classical mechanics appears as a strange intermediate case. 
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5 Conclusions 



The identities (pCi|), (p3l|), ( p2| ) present the main result of this paper. The show that the matrix elements 
of arbitrary operator in quantum Toda chain can be expressed with help of finitely many integrals which 
possess remarkable properties. 

There is a difference with what we have in Sine-Gordon theory in infinite volume. Indeed, the matrix 
elements for Toda chain are given by integrals of arbitrary deformed differentials with respect to fixed half- 
basis of deformed cycles. In Sine-Gordon case we took arbitrary half-basis. The reason for that is the 
difference in the type of reality conditions. 

In this connection it would be very important to consider the deformation of Toda chain with trigono- 
metric R-matrix and more complicated reality conditions which is more close to Sine-Gordon case. We hope 
that in this situation there will be a complete duality between deformed differentials and deformed cycles. 
Acknowledgments. I am grateful to O. Babelon, D. Bernard, L. Faddeev, E. Frenkel and N. Reshetikhin 
for discussions. I would like to thank RIMS at Kyoto University where the work was finished for hospitality. 
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